Semiclassical instability of dynamical warp drives 



Stcfano Finazzij^'B Stcfano Libcrati.^'0 and Carlos Barcel6^'[|| 

^SISSA, via Beirut 2-4, 34151 Trieste, Italy and INFN sezione di Trieste, via Valeria 2, 34127 Trieste, Italy. 
^Instituto de Astrofisica de Andalucia, CSIC, Camino Bajo de Huetor 50, 18008 Granada, Spain 

(Dated: July 14, 2009) 

Warp drives are very interesting configurations in general relativity: At least theoretically, they 
provide a way to travel at superluminal speeds, albeit at the cost of requiring exotic matter to exist 
as solutions of Einstein's equations. However, even if one succeeded in providing the necessary exotic 
matter to build them, it would still be necessary to check whether they would survive to the switching 
on of quantum effects. Semiclassical corrections to warp-drive geometries have been analyzed only 
for eternal warp-drive bubbles traveling at fixed superluminal speeds. Here, we investigate the more 
realistic case in which a superluminal warp drive is created out of an initially flat spacetime. First of 
all we analyze the causal structure of eternal and dynamical warp-drive spacetimes. Then we pass to 
the analysis of the renormalized stress-energy tensor (RSET) of a quantum field in these geometries. 
While the behavior of the RSET in these geometries has close similarities to that in the geometries 
associated with gravitational collapse, it shows dramatic differences too. On one side, an observer 
located at the center of a superluminal warp-drive bubble would generically experience a thermal 
flux of Hawking particles. On the other side, such Hawking fiux will be generically extremely high 
if the exotic matter supporting the warp drive has its origin in a quantum field satisfying some form 
of quantum inequalities. Most of all, we find that the RSET will exponentially grow in time close 
to, and on, the front wall of the superluminal bubble. Consequently, one is led to conclude that the 
warp-drive geometries are unstable against semiclassical backreaction. 

PACS numbers: 04.20.Gz, 04.62. +v, 04.70.Dy 
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I. INTRODUCTION 

Since Alcubierre introduced them [l[ , warp drives have 
been certainly one of the most studied spacetime geome- 
tries among those requiring exotic matter for its existence 
(see for a recent review about "exotic spacetimes"). 
They were immediately recognized as interesting configu- 
rations for two main reasons. First, they provide, at least 
theoretically, a way to travel at superluminal speeds: A 
warp drive can be described as a spheroidal bubble sep- 
arating an almost flat internal region from an external 
asymptotically flat spacetime with respect to which the 
bubble moves at an arbitrary speed. The idea is that 
although on top of a spacetime nothing can move with 
speeds greater than that of light, spacetime itself has no 
a priori restriction on the speed with which it can be 
stretched. Second, they provide an exciting ground for 
testing our comprehension of general relativity (GR) and 
quantum field theory in curved spacetimes (for instance 
when investigating warp-drive implications for causality 

. . 

In its original form [1| the warp-drive geometries are 
described by the simple expression 

ds^ = -c^dt^ + [dx - v{r)dtf + dy^ + dz'^ , (1) 
where r = ^ [x — vot)^ + + is the distance from 
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the center of the bubble and vq the warp-drive velocity. 
Here and thereafter v = vofir) with / a suitable smooth 
fimction satisfying /(O) = 1 and /(r) — > for r oo. 

After the warp-drive spacetimes were proposed, the 
most investigated aspect of them has been the amount 
of exotic matter (i.e. energy-conditions-violating matter) 
that would be required to support them.^ It was soon 
realized that this was not only related to the size of the 
warp-drive bubble but also determined by the thickness 
of the bubble walls 0]. It was found that, if the ex- 
oticity was provided by the quantum nature of a field, 
satisfying therefore the so-called quantum inequalities 
(QI),^ then the violations of the energy conditions would 
have to be confined to Planck-size regions, making the 
bubble-wall thickness A to be, accordingly, of Planck 
size [A < 10^ {vq/c) Lp, where Lp is the Planck length]. 
However, it can be shown that very thin walls require 
very large amounts of exotic matter: e.g., in order to 
support a warp-drive bubble with a size of about 100 m 
and propagating at vq ~ c, one would need a total nega- 
tive energy \E\ S> IO^^A/q.^ Perspectives for warp-drive 



^ While initially it was supposed that exotic matter was needed 
only for superluminal warp drives (vo > c), it was later rec- 
ognized Q, [3I that energy-conditions- violating matter is needed 
also for subluminal speeds. This points out that the need of 
exotic matter is peculiar of the warp-drive geometry itself, not 
appearing only in the superluminal regime. 

^ See 01 for a review about quantum inequalities applied to some 
exotic spacetimes. 

^ If one could somehow avoid the QI, then it would be possible to 
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engineering can be improved by resorting to a modified 
warp-drive configuration with a reduced surface area but 
the same bubble volume The total amount of neg- 
ative energy required to support these warp drives be- 
comes quite small (for example, \E- \ « 0.3Mq for a 100 
m-radius bubble, although one has to add as well some 
positive energy outside the bubble £'+ « 2.5Af0), bring- 
ing the warp drive closer to a realistic solution albeit still 
far from foreseeable realizations. 

Regarding the feasibility of warp-drive configurations, 
a parallel line of research has focused on the study of their 
robustness against the introduction of quantum correc- 
tions to GR. In particular, in Q it was studied what 
would be the effect of having semiclassical corrections in 
the case of an eternal superluminal warp drive. There, 
it was noticed that to an observer within the warp-drivc 
bubble, the backward and forward walls (along the di- 
rection of motion) look, respectively, like the horizon of 
a black hole and of a white hole. By imposing over the 
spacetime a quantum state which is vacuum at the null 
infinities (i.e., what one may call the analog of the Boul- 
ware state for an eternal black hole) it was found that 
the renormalized stress-energy tensor (RSET) diverges 
at the horizons.^ Independently of the availability of ex- 
otic matter to build the warp drive in the first place, the 
existence of a divergence of the RSET at the horizons 
would be telling us that it is not possible to create a 
warp-drive geometry within the context of semiclassical 
GR: Semiclassical effects would destroy any superlumi- 
nal warp drive. However, in a more realistic situation, 
a warp drive would have to be created at a very low ve- 
locity in a given reference frame and then accelerated 
to superluminal speeds. One may then expect that the 
quantum state globally defined on such dynamical geom- 
etry would be automatically selected by the dynamics 
once suitable boundary conditions are provided (e.g., at 
early times). This is indeed the case for a gravitational 
collapse where it can be shown that, whenever a trap- 
ping horizon forms, the globally defined quantum state 
that is vacuum on has to be thermal at (at the 
Hawking temperature) and regular at the horizon. In 
other words, the dynamics of the collapse avoids select- 
ing a Boulware-like state, with its associated divergence 
at the horizon, ending up instead selecting the analo- 
gous for collapsing configurations of the Unruh vacuum 
state defined on eternal black holes, which leads to a per- 
fectly regular RSET. Is the dynamics of the creation of 
a warp drive, with its associated selection of the global 
vacuum state, able to avoid the presence of divergences in 
the RSET? Indeed, in 0] it was already noticed that an 
Unruh-like state rather than a Boulware-like state should 
be expected to describe the quantum state characterizing 



built warp drives with much larger wall thickness. For example, 
for A ~ 1 m, one would only need \E\ > 1/4 Mq. 
^ Throughout this paper we shall work in the Heisenberg represen- 
tation so that only operators, not the states, evolve in time. 



a superluminal warp-drive creation. 

In this paper we want to settle this issue by explicitly 
considering the case of a warp drive which is created with 
zero velocity at early times and then accelerated up to 
some superluminal speed in a finite amount of time. This 
can be viewed as the warp-drive analog of a semiclassical 
black hole collapse [l3|. By restricting attention to warp 
drives in 1 -f 1 dimensions (since this is the only case for 
which one can carry out a complete analytic treatment), 
we have calculated the RSET in the warp-drive bubble. 
As we were expecting, we find that in the center of the 
bubble there is a thermal flux of particles at the Hawk- 
ing temperature corresponding to the surface gravity of 
the black horizon. However, the surface gravity can be 
shown to be inversely proportional to the thickness of the 
bubble walls which, as said, has to be of the order of the 
Planck length if some form of QI holds. Hence, one has 
to conclude that an internal observer would soon find it- 
self in the uncomfortable condition of being swamped by 
a thermal fiux at the Planck temperature. Even worse, 
we do show that the RSET does increase exponentially 
with time on the white horizon and close to it (while 
it is regular and small on the black one). This clearly 
implies that warp drives would become rapidly unstable 
once superluminal speeds are reached. 

The plan of the paper is the following. In Sec. [H] we 
study the causal structure of both an eternal and a dy- 
namic warp drive. In Sec. Illll we discuss the propagation 
of light rays in a dynamical warp-drive geometry in close 
analogy with the standard treatment for black holes. ^ 
Finally, in Sec. [IV]we calculate the full RSET using the 
technique adopted in [l0| for black hole formation and 
look for its divergences. A summary of our results is 
given in Sec. fVl 



II. CAUSAL STRUCTURE OF A 
SUPERLUMINAL WARP DRIVE 

We investigate the causal structure of a warp drive, 
following the method presented in [T^ . for spacetimes 
whose metrics can be written in Painleve-GuUstrand co- 
ordinates. We begin from an eternal warp drive, moving 
at constant velocity, and then we study a dynamic situ- 
ation in which a warp-drive bubble is accelerated. 



A. Eternal superluminal warp drive 

In 1+1 dimensions Alcubierrc's metric Eq. (jT]) reduces 

to 

ds'^ = -c^dt"^ + [dx ~ v{r)dtf . (2) 



A similar result was found by Gonzalcz-Diaz 
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Here r is defined as the signed distance from the center 
of the bubble, r = x — VqI. Again we define v(r) = Vof^r) 
but / is now taken to be defined also for negative values of 
r. Its boundary conditions will be /(O) = 1 and /(r) — > 
for r — > ±cx). For illustrative purposes, let us choose the 
following simple bell-shaped function: 

cosh [r/aj 

and vq > c, that is, our warp drive is superluminal. We 
want to stress that all the results of this paper do not 
depends on the particular choice of the function /. They 
are still valid providing that it satisfies the above condi- 
tions. We have chosen a particular form for it just for 
simplicity. Now, by using [t, r) coordinates the metric 
Eq. ^ reads 
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FIG. 1: Lines of constant u (solid lines) and w (dashed lines) 
for the Alcubierre warp drive, a = c= l,Q = 2. 



-c^df + [dr - v{r)dt] 



v{r) = v{r) ~ Wo . 

(4) 

Note that v < because the warp drive is right going 

{v > 0) but v{r) < Vq. 

By definition, at any time the center of the bubble is 
located at r = 0. The causal structure of this metric can 
be analyzed following [l^ . Let us define the warp-drive 
Mach number a = vq /c. The shift velocity becomes 



v[r) 



ac 



cosh (r/a) 



- 1 



(5) 



Two horizons appear when a > 1. Their positions are 
found by putting v equal to — c, 



ri,2 



Taln(/3+ 



/3: 



a - 1 



>1. (6) 



Because of our simple profile choice, the two horizons are 
symmetrically located with respect to r = 0. In more 
general situations ri and r2 will be completely arbitrary 
satisfying only ri < r2. 

Right- and left-going null coordinates u and v can be 
defined as 



du = dt " 

dw = dt + 



dr 



c + v{r) ' 

dr 
c — v{r) 



(7) 
(8) 



We note that the spacetime is divided into three distinct 
regions: I (r < n), II (n < r < r2), and III (r > r2). 
There are no u rays connecting these regions, while w rays 
cross all the regions; when r approaches the horizons, u 
diverges logarithmically. Integrating these equations we 



find 



ui = t 



c(a- 1) 
a/3 



r 

Ull=t + 



:ln 



-{r-ri)/a _ ^ 



c(q! - 1) 

a/3 



c(a- 1)7^231 
r 

Mm = t + 



In 



-{r-r2)/a _ I 



-(] — ri)/a 



c (a — 1) 
a/3 



c(a - 1) ^/32 - 1 
w = t + 



:ln 



-(r-r2)/a _ I 



-(] — ri)/a 



1 — Q-{r-r2)/a 



r 



c(a + 1) 
2a7 



: arctan 



7^ 



(9) 



(10) 



(11) 



(12) 



where /3 is defined in Eq. ([6]) and 7 = a/ (a + 1) < 1. In 
Fig. [1] we plot the lines of constant u and w. The horizon 
at r = ri corresponds to ui — > -l-oo and to uu —^ +00, 
while the horizon at r = 7'2 corresponds to uu —> —00 
and to uni — —00. 

We define a signed surface gravity on our two horizons: 



_ dv{r) 



dr 



^ i^i^^^il^^ . ±. . (13) 



a/3 



So doing, the surface gravity associated with the first 
horizon is positive ki = k > 0, while the one associated 
with the second horizon is negative K2 = < 0.^ As we 



In the definition of surface gravity, more correctly one siiould 
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noticed for the positions of ri and r2, in a general situa- 
tion in which / is not symmetric, the two surface gravities 
may have different absolute values. However, the one as- 
sociated with the first horizon (respectively, second hori- 
zon) will be always positive (respectively, negative). We 
will soon show that these two horizons represent a black 
and a white horizon, respectively. Here on, we will con- 
sider the two surface gravities to have the same absolute 
value. Retaining different absolute values will not lead 
to more general results, but will just make the notation 
heavier. 

As in Sec. 6 of [l^l, one can find some transformation 
Ui = Ui{ui) that close to the horizons behaves as 

Ui{ui ^ +00) ~ C/i + Aie-™' , (14) 
Un{un ^ ±^) ^ Un± T Anie^'^"" , (15) 
Um{uui ^ -00) ^ Um - Ame™'" , (16) 

where [/j, C/ii_|_, J7ii_, and [/m are arbitrary constants 
and Ai, An+, An_, and Am are positive constants. It is 
possible to choose these transformations such that they 
match on the horizons, obtaining a global U varying reg- 
ularly from 4-00 to — cx) as r varies from ~cx) to -f-oo. This 
matching will become natural when dealing with dynam- 
ical configurations in which the warp drive is created by 
accelerating the bubble from an initial zero velocity. In 
order to do that we can choose Ui = Uu^ = Ubh and 
Uu- = Uiu = ?7wH < C^BH- In this way the three trans- 
formation have the following domains: 

uiG(-oo,+oo) [/i e (+CX), C/bh) , (17) 

uii e (+00, --00) Uu G {Ubh^Uwh) , (18) 

Will e (-00,-1-00) =^ f/iii G (t/wH, -00) . (19) 

For instance, specific transformations having the re- 
quired properties are the following: 

Ui^l+ e-^"' , (20) 

C/ii = itanh(^) , (21) 

f/iii = -^-e'^""V (22) 

Now we can bring the right and left infinities to a finite 
distance by using a compactifying transformation like 

Ui = arctan(t/i) , (23) 
Um = arctan(;7iii) , (24) 
W = arctan(w) . (25) 

The Penrose diagram for this spacetime is plotted in 
Fig. [51 Notice that the diagram does not correspond to a 




FIG. 2: Penrose diagram of an eternal warp drive. Lines 
of constant r (solid lines) and of constant t (dashed lines). 
Future and past horizons a,t r — ri,r2 (heavy dashed lines). 
The geometry can be extended to the future of Jff^^ and J^^l 
and to the past of Jif(^ and Jif£~ . 

maximal analytical extension but to a particular patch of 
the total spacetime. The dashed lines signal the locations 
at which the geometry can be extended. The two exter- 
nal regions r > r2 and r < ri appear to an observer living 
inside the bubble as an eternal white hole and an eternal 
black hole, respectively. For this reason, we call and 
respectively, the black horizon and the white hori- 
zon of the warp drive. Notice however that for the inner 
observers in an eternal configuration, at ri there is also 
a white horizon [J^^(t — — oo,r = ri) in the diagram] 
and at 7*2 a black horizon [J/iCq [t = -|-oo,r = r2) in the 
diagram] . The geometry can be extended through these 
two null lines. We do not picture the extended regions 
as there arc ambiguities in the prescription of the matter 
distribution in those other "universes."'' 



B. Dynamic superluminal warpdrive 

What happens to the causal structure when we con- 
sider the creation of a superluminal warp drive starting 
from initially flat spacetime? For concreteness, we study 
a simple case in which we reach the final velocity v = vq 
at a finite time which we take to be t = 0. We modify 



multiply the velocity derivative by c to get an acceleration. How- 
ever we use this slightly modified definition to avoid the appear- 
ance of too many constant c factors in our formulas. 



In the analysis of the maximal extention of the Schwarzschild or 
Kerr spacetimes this problem is not present as one considers only 
vacuum solutions of the Einstein equations. 
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the metric Eq. Q introducing a switching factor S(t): 

ds^ = -c^dt^ + [dx - w(r, t)dtf , (26) 

where 

v{r,t) = vo5{t)f{r) , (27) 
with f{r) defined in Eq. ([3]) and 

'' e*/^ if t < , 



1 



if t>Q . 



(28) 



Again we can change coordinates, keeping the center of 
the bubble at rest (r = 0). This can be obtained by 
defining 



dr = dx ~ VQd{t)dt 



(29) 



This is an exact differential form and can be integrated 
to get 




1 



if t<0 , 
if t>0 . 



(30) 



Again, the metric becomes 



ds'^ = -c^dt'^ + [dr - v{r,t)dtf 
vir,t)^voSit) [f{r)-l] 



(31) 



and the apparent horizons associated with the t slicing 
are found by setting v — — c. In this case a solution does 
not exist for any value of i, so that the apparent hori- 
zons are created at infinity at some finite in. We show 
this below. Let us write the equation for the apparent 
horizons in the following form: 



/M = 1 - 



(32) 



The function / takes all the values between and 1. In 
particular, /(r) for r — > ±oo and /(O) = 1. Besides, 
the right-hand side of Eq. \i2\ is a monotonic function 
of t. such that, for t —oo, 1 — c/ {vo5{t)) — oo 
and reaches the value 1 — c/vq > for t > 0. As a 
consequence, there exists a time tn < Oso that for t > tn 
there are always two solutions of Eq. (j32p , corresponding 
to a black and a white horizon. These horizons are born 
aX t ^ tu at spatial infinity and at t = they settle at 
two fixed positions ri and r2. 

Keeping these points in mind, we are able to build the 
Penrose diagram for the dynamic warp drive (Fig. [3]). 
At early times the metric is approximately Minkowskian, 
because 5{t) — > for t —oo. Therefore, the causal 
structure is initially Minkowskian. Then, it progressively 
changes till t = 0. At this time one has built a stationary 
warp drive, just as in the previous section. After this 
time, the Penrose diagram looks exactly equal to that in 
Fig. [21 The final Penrose diagram is just obtained by 




FIG. 3: Penrose diagram of a dynamic warp drive. Lines 
of constant r (solid lines) and of constant t (dashed lines). 
The lines of constant r become null at the apparent horizons 
(heavy dashed lines). 



gluing together the two behaviors. Again, we are not 
drawing an analytically extended diagram but only the 
relevant patch for the analysis that follow in this paper. 
Given that a timelike observer can reach J^^^ and 
in a finite proper time, the geometry can be extended in 
the future, beyond these lines. This is actually a sub- 
tle point. In fact, and (which are linked to 
the formation of a white horizon) are on the boundary 
of the Cauchy development of J"' . In this sense, they 
are Cauchy horizons given that initial data are assigned 
only on J'^ . Hence, as noticed in Sec. Ill Al an eventual 
extension would not be unique. In any case such an ex- 
tension will not be relevant for what will follow, given 
that we shall limit ourselves to investigating the asymp- 
totic behavior of the RSET associated with the onset of 
the super luminal warp drive. 

Let us highlight here that the dynamic way we have 
used here to create the warp drive is not the only possible 
way. In Appendix \^ we will also use a different interpo- 
lation between Minkowski and the warp drive in which 
the horizons appear also at finite time but at finite r 
positions (similar to what happens when a homogeneous 
star collapses to form a black hole). 



III. LIGHT-RAY PROPAGATION 

The just discussed causal structure of the dynamical 
warp drive is naturally leading to the expectation that 
some sort of Hawking radiation will be produced in a su- 
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pcrluminal warp drive (as well as some transient particle 
emission). It is well known (see, for example, [13|) that 
all the information about particle production is encoded 
in the way in which light rays propagate in a spacetime. 
That is, it is enough to know how light rays are bended 
in order to analyze the phenomenon of particle creation. 
In the dynamical warp drive there is a single past null 
coordinate but three different future null coordinates as- 
sociated with the final regions I, II and III, as described 
before. From here on, we will be dealing exclusively with 
the connection between the past null coordinate U at 
and the future null coordinate uu at J^q in the interior 
of the bubble. Therefore, we will use u to denote uu 
whenever this does not lead to confusion. As discussed 
in the relation U = p{u) encodes all the relevant 
information about particle production. 

We want to study the features due to the two main 
properties of a dynamical warp-drivc geometry, i.e., the 
spacetime is Minkowskian at early times and it is a warp 
drive at late times. In particular, we are not interested 
in the transient features depending on how the transi- 
tion between these two regimes is performed. Namely, we 
need only the behavior close to the horizons and at late 
times inside the whole bubble. It is clear from Fig.[3]that, 
if one stays at constant r inside the bubble and moves 
forward in time one crosses u rays which pass closer and 
closer to the black horizon. Therefore, once we have de- 
termined the behavior of p{u) close to the horizons,^ we 
automatically also have the required information at late 
times in the whole bubble. 

In general, the relation U ~ p{u) is obtained by inte- 
grating the differential equation for the propagation of 
right-going light rays 



dr 
'dt 



v{r,t) . 



(33) 



Note that while in the previous section we considered a 
specific form of v{r,t) [Eq. (PT|) ] in order to discuss the 
causal structure of the associated spacetime, here (and in 
what will follow) the discussion will hold for any v{r,t) 
that satisfies the requirements v{r,t) — > 0, for t — > — oo 
(sufficiently rapid for the spacetime to be asymptotically 
flat), and v{r,t) = v after some finite time and within 
the warp-drive bubble. Given the assumption that the 
spacetime settles down to a stationary warp-drive con- 
flguration, at late times the velocity profile will depend 
only on the r coordinate. We can write, as in the sta- 
tionary case Eq. ([7]) 



du = dt 



dr 



- v{r) 



(34) 



To find the required asymptotic relation one has to in- 
tegrate this equation in the limit r — > ri^2- There, the 



Appendix |A] gives a specific example for which this relation can 
be computed exactly in the whole spacetime. 



velocity can be expanded as 

w = -c±K(r-ri,2) + C»((r-ri,2)^) . (35) 



Thus, we obtain 



u ~ t ^ — In |r — ri 2I 



(36) 



On the other hand, the coordinate [/, obtained by inte- 
grating Eq. ([33|) at early times, reduces to the Minkowski 
null coordinate 



U{t ~oo) 



(37) 



and is regular in the whole spacetime, in particular, on 
the horizons. For instance, on a fixed t slice in the sta- 
tionary region, we can write [/ as a regular function of 
r 



f/± =W±(r-ri,2) 



(38) 



where wc denoted with ?7_|_ (respectively, t/_) the specific 
form of U close to the black (respectively, white) horizon, 
and U± are analytic functions. Inserting Eq. ([36)1 in the 
above expression, at the same fixed time, 



U± = p{u ±00) = 7'±(eT'"") , 



(39) 



where V± are again analytic functions. Note that the 
forms of these functions do not depend on the particu- 
lar time slice chosen to perform the matching between 
Eqs. (|36p and (|38p . In the proximity of the stationary 
horizons u — > ±00, so e^''" — ^ and the function p can 
be expanded around the horizons. Up to the first order 
we get 

U ^plu^ ±00) = C/bh T ^±e=F'"" + O (e=F2«") , (40) 

WH 

where A± are positive constants. 

This is indeed the asymptotic behavior one would ex- 
pect in the presence of trapping horizons and, for the 
black hole case, it is the standard relation between u and 
U . In fact, it leads to the conclusion that an observer 
at will detect Hawking radiation with temperature 
njl-n. It is important to note that the result is completely 
general. The asymptotic behavior of f/ = p{u) for large 
absolute values of m, which is the only feature of p{u) rel- 
evant for the analysis of this paper, does not depend on 
the specific velocity profile adopted. It is only necessary 
that it interpolates from Minkowski spacetime at early 
times to a stationary warp-drive geometry at late times. 

While Eq. (|40|) is exactly of the expected form, its im- 
plications for particle production arc not as straightfor- 
ward as in the black hole case. In fact, in the warp-drive 
geometry the late-time modes labeled by u will not be 
standard plane waves in an asymptotically flat region of 
spacetime as they will be characterized by the strange 
form given in Eq. (|36p . Of course, if the surface gravity 
K is large enough so that the typical wavelength of the 
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emitted radiation is much smaller than the bubble size, 
then the plane- wave approximation is fine and in the cen- 
ter of the bubble, at late times, one will measure standard 
Hawking radiation at temperature T. Nonetheless, in the 
general case, even if the calculation for the Bogoliubov 
coefficient is the standard one the particles created 
will not be standard plane waves. The physics associ- 
ated with the particle production by the white horizon 
is even less clear. This is why, in the next section, we 
shall consider the behavior RSET to get more significant 
information. 

In order to do so, we shall also need the relation be- 
tween ingoing and outgoing left-going rays. In fact, these 
modes are excited too when the warp drive forms, even 
if we do not have a thermal particle production as for 
right-going modes. Left-going rays are the solution of 
the following differential equation [see Eq. ([33]) ]: 



dv 

- = -c + v{r,t). 



(41) 



Looking at Fig. [31 we note that left-going rays do not see 
the horizons, that is, they cross them from to . 
As a consequence, both the past and the future null co- 
ordinates W and w are defined at the asymptotic region 
outside the bubble. However, after the geometry inside 
the bubble has settled down to its final stationary form, 
we can define, just for convenience, a different coordinate 
w inside the bubble, as in Eq. ([8]): 



dw = dt 



dr 



c — v{r) 



(42) 




Note that w and w may or may not coincide depending on 
how fast the metric in the external region settles down to 
its final stationary form (refer to Fig. |4] and Appendix [XI 
for an example in which they do not coincide). 

W is obtained in the usual way, by integrating Eq. (|4T|) 
at early times, when the spacetime is Minkowski 



W{t -oo) 



r 

t + - 
c 



(43) 



The relation W = q{w) can be found explicitly for specific 
cases (see Appendix [A|. The important point is that one 
can prove that this relation is always regular, so that 
it cannot give place to any phenomenon like Hawking 
radiation. Of course, one can choose a very unusual way 
to interpolate from Minkowski to the warp drive, such 
that a lot of particles are created in this sector, but this 
is not a general feature of dynamical warp drives. If we 
use a regular enough transition, only transient effects are 
present in this sector. In Appendix [Xl we shall show that 
it is possible to find such a transition. 

As an example, in Fig. [31 we plot both right-going and 
left-going light rays propagating in the particular dy- 
namic warp-drive spacetime analyzed in Appendix [^ 



FIG. 4: Light rays propagating rightward {solid lines) and 
leftward (dashed lines) in the plane {t, r) in a warp-drive 
spacetime with velocity profile of Eq. (|A2|l . The out re- 
gion in which the geometry is a stationary warpdrive is at 
r < ±arccosh(t-|- 1) (heavy solid lines). At t < the metric is 
Minkowskian. The horizons at ri and r2 (heavy dashed lines) 
are formed at Th = 1. Please refer to Appendix[Xlfor details. 



IV. RENORMALIZED STRESS-ENERGY 
TENSOR 

For the calculation of the RSET inside the warp-drive 
bubble we can use the method used in jT3| for a collapsing 
configuration to form a black hole. In null coordinates U 
and W, affine on J^~, the metric can be written as 



ds^ 



'C{U, W)dUdW . 



(44) 



As we described in the previous section, in the out re- 
gion (the region in which the metric is stationary, i.e., 
V depends only on r) we can also use a different set of 
null coordinates m, w. The coordinate u is affine on M'f^ 
and w is the coordinate defined in Eq. In these 
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coordinates the metric is expressed as 
ds^ ~ —C{u,'w)dudw , 

which imphes 

C{u, w) 



C{U, W) = 



p{u)q{w) 



and 



U=p{u), W^q{w) 



(45) 



(46) 



(47) 



These transformations are such that the C has precisely 
the form of the future stationary warp-drive geometry, 
that is, it depends only on r through 

For concrctencss let us refer to the RSET associated 
with having a single quantum massless scalar field living 
on the spacetime. In this case the RSET components 
have the following form : 



iZTT 

Tuw = Twu = T^CR . 



Tuu 
Tww 



(48) 
(49) 
(50) 



Qualitatively, things would not change if there were other 
fields present in the theory. The only modification will 
be that the previous expressions will get multiplied by 
a specific numerical factor. It is clear that in the in re- 
gion (where the spacetime is Minkowskian) the RSET is 
trivially zero. 

Now, for a metric in the form of Eq. (jH]), the curvature 
can be calculated as in [l^: 



i? = DlnlCI = {-gy^'^d^ [(-g)i/V"a,ln|C| 

4 

Replacing this result in Eq. ([50|) we obtain 

1 



'-WU 



247r 



dudw\^\C\ 



Using transformations pG]) and (j47|) in Eqs. (|48 
and (1521) we obtain 



Cl/292 C-1/2 _ gl/2a2 .-1/2 
1 1 



(51) 
(52) 



Tww 



1 1 r 
12^^ 
1 1 
12^^ 

Tuw = Twu = — 



247r -pq 



dud^ In |C| 



(53) 
(54) 
(55) 



We can express the derivatives with respect to u and 
w in terms of derivatives with respect to r and t. Note 
that the following expressions are obtained from Eqs. ([M]) 
and (I42p , so they arc valid only when the velocity profile 



depends only on the r coordinate. In this fashion, we can 
study the stress-energy tensor at the end of the creation 
of the warp drive (or at the end of a collapse if we are 
studying a black hole). 



Ur Wr 

ut lit 



-l/{c + v) l/{c-v) 

1 1 



(56) 



Inverting the derivative matrix we obtain the required 
result: 



du 



2c 



-dr 



2c 



-a. 



c + V 

C — V 



dt 



2c 



dt 



(57) 
(58) 



Please refer to Appendix [B] for more details. 

We are interested in calculating the RSET inside the 
bubble when the two horizons have been formed and 
when the configuration has settled down to a station- 
ary warp drive. Since in this region the velocity depends 
only on r, we can replace the derivatives in Eqs. (|53|1. 
(I5il) and ((53 with 



du 
d,r, 



dr 



c 



dr 



(59) 
(60) 



Here, we have put c = 1, for simplicity, and we have used 
Eq. (HSI, C = 1-v'^. Moreover we indicate with ' the 
differentiation with respect to r. After some calculations 
we obtain 



(51/2^2(^-1/2 ^ ^1/2^2^-1/2 ^ _^ 



{C'f -2CC" 



pV2a2^-V2^_L[3p2-2pp] 

e^'dlq-^/^ = ^[3f-2q-ci-] 



duduM\C\ 

Using again C = 1 — we get the final result: 



(61) 
(62) 

(63) 

(64) 



T 



uu 



WW 



f-uw 



1 1 

'487r^ 

1 1 

'48^^ 



«'2 + (l-w2) vv" 



Z;'2 + (1~Z;2) vv" 



3p2 - 2p p 



3g2 - 2q q 



Twu = ^75-— (1 
AViTT pq ^ 



(65) 

(66) 
(67) 



One can check that these quantities do not diverge at the 
horizons, just like in [l3|. However we want to look at the 
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energy density inside the bubble and try to understand measured by a set of free-falling observers, whose four 

whether it remains small or not as time increases. velocity is simply = (l,w) in {t,r) components. 
In particular, it is interesting to look at the energy 

I 



These observers measure an energy density p: 



L^^u,^ u,^ ^ in-t 2vTtr + v'^Trr — U^Tuu + 2UtWtTuw + W^Tww 



T,,„u'iu'i ^ Tff 



+ 2v [UtUrTuU + [UtWr + WtUr) Tuw + WtWrTww] + v"^ [UrTuu + 2UrWrTuW + W^Tww] 
= {Ut + vUrf Tuu + 2 ([/t + vUr) (Wt + vWr) Tuw + [Wt + vWrf Tww 



(l + v) 1-w 



:Tuw 



(1 



tTww 



1 

48^ 



'2{v^ 



r.' 2 



4:V 



f{u) 



(1 



(1 



(68) 



where we have defined 



giw) 



w) - 2q{w) 'q{w) 



(69) 
(70) 



We want to study what happens to a spaceship placed 
at rest in the center of the bubble to investigate whether a 
warp drive can be used as a transportation device. More- 
over we want to check whether the components of the 
RSET in a regular coordinate system are regular at the 
horizons. Freely falling observers are at rest in the center 
of the bubble but left moving otherwise. 

Looking at the above expressions for any component 
of the RSET or at the energy density p, we sec that they 
can be split as a sum of three terms, one purely static, 
depending only on the r coordinate through the shift 
velocity v, and two dynamic pieces depending also on 
the u and w coordinates, respectively. They correspond 
to energy traveFling on right-going and left-going light 
rays, respectively, eventually red/blue-shifted by a term 
depending on r. Being interested in studying the energy 
density measured by free-falling observers, we write 



P — Pst + Pdyn-u + Pdyn- 



(71) 



where we defined static terms (labeled by subscript st) 
and dynamic terms (labeled by subscript dyn) for each 
component. Right-going and left-going terms are respec- 
tively labeled by u and w. 



Pst 



1 

24^ 



-I 2 
-V 



2v 



Pdyn—u 
Pdyn- 



_ 1 fiu) 



487r + 

1 gjw) 

487r (1 - v) 



,2 ' 



2 ■ 



, (72) 
(73) 
(74) 



Let us start with pdyn-is ■ Its denominator is bounded for 
each r because the shift velocity is negative, and its nu- 
merator is vanishing with time. It is easy to convince one- 
self that all the contributions to the RSET coming from 
the w sector are not universal but depend exclusively on 
the specific interpolation between the early Minkowski 
spacetime and the final warp-drive spacetime. It is al- 
ways possible to choose an interpolation so that all these 
contributions vanish at late times (see Appendix E]) . 

From now on, we neglect the dynamic term pdyn-w 
and study 



P — Pst + Pdyn-u 



(75) 



A. RSET at the center of the warp-drive bubble 

We shall now study the behavior of the RSET in the 
center of the bubble at late times. In this point v(r = 
0) — v'{r ^ {)) = and the static term pst vanishes. 
By integrating Eq. (j34[) in the stationary region, one can 
show that u{t,r) is linear in t so that, for fixed r, it 
will acquire arbitrarily large positive values (see Fig. [3]). 
The dynamic term in Eq. (|73p can be evaluated by using 
a late-time expansion oi p{u) [see Eqs. (|40p and (|C2[) ]. 
then plugged in Eq. This is derived in Appendix [Cl 

and it is equal to 



A 



2+ 



A, 



-2ku 



+0(e-3'"")} . (76) 



With this expansion it is easy to see that for u 

P ■ 



K 

487r 



(77) 



This result may be easily understood in the following way. 
The surface gravity of the black horizons is the veloc- 
ity derivative evaluated at this horizon k = {dv / dr)r=r-i- 
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Moreover, the energy density of a scalar field at some 
finite temperature T in 1 + 1 dimensions is simply 



PT 



UJ dk rp2 

Ty2^ ^ 12 



(78) 



Defining the Hawking temperature in the usual way as 
Th = k/2ti, it is easy to see that Eq. ((77|) can be rewrit- 
ten exactly as p = {t: /12)Th^ . Thus, the observer inside 
the warp-drive bubble will indeed observe thermal radi- 
ation at the temperature Th- 



B. RSET at the bubble horizons 

Let us move now to study p close to the horizons. Note 
that both pst and pdyn-u are divergent at the horizons 
(r — > Ti 2) because of the (1 -f v) factors in the denom- 
inator. Just like in Oil, for a black hole, these diver- 



gences exactly cancel each other, but something different 
happens at the black and white horizon. Expanding the 
velocity up to the second order one gets 

v±{r) = -1±K (r - ri,2) + ^cr [r - 2)^+0 {{r - ri,2)^) 

(79) 

where cr is a constant. 

Close to the horizons the static term of Eq. ([7^ then 
becomes 



Pst- {t ~ ri 2) = 



(r-ri,2) K(r-ri,2) 



-0{l). 
(80) 

Similarly one can expand the dynamic term pdyn-u to 
the same order. This involves determining the function 
,f{u) in the proximity of the horizons. Indeed it is not 
difficult to show that (see Appendix [C]) 



lim fiu) = ll 

r — >r-i- I 



2± 



'A 



(r-ri,2)' + 0((r-ri,2)') 



(81) 



Inserting this form in Eq. ((73|) we obtain: 



C. RSET at late times approaching M'fj 



Pdyn-u [r ~ ri,2) 



487r {r-ri2f^'^('r-ri.2) ^^^^^ 

(82) 

It is now clear that the total p is 0(1) on the horizon and 
does not diver ge ( as expected from the Fulling-Sweeny- 
Wald theorem [l5|). 

However, let us look to the subleading terms 



In Sec. HVB] we studied the RSET on the black and 
white horizons. However, it is interesting to analyze bet- 
ter its behavior as it approaches the Cauchy horizon . 
As one can see from Fig. HJ every u ray reaches values of 
r very close to r2, at sufficiently late times. This means 
that, even for large positive values of u, some time exists 
after which the approximate behavior of all the u rays is 
well described by 



p{r 



'''■'^ ^ 48^ 



A 



2± 



- 2 



A3± 



Ai±J Ai± 

+ C±+0(r-ri,2) , (83) 



where C± are constants. We can easily see that the be- 
havior close to the black horizon is completely different 
from that close to the white horizon. In the former case 
the energy density as seen by a free-falling observer is 
damped exponentially with time. In the white horizon, 
however, this energy density grows exponentially with 
time. This means that, moving along , p is large and 
negative diverging while approaching the crossing point 
between and J^q . 

This asymptotic divergence is physical and not a mat- 
ter of selection of coordinates. In a very short time after 
the white horizon is formed (of the order of 1/k), the 
backreaction of the RSET in this region of spacctinie is 
no longer negligible but rather very strong. 



so that 



t + -\n (7-2 



r2 



r) , 



(84) 



(85) 



Keeping in mind this point, let us study the term pdyn-wj 
at fixed u, when time increases. There is a time at which 
the u ray will be close enough to r2, such that the de- 
nominator in Eq. (j73p can be approximated by 



(86) 



(l + vy c::^ K^ir2-r)- , 
which becomes, thanks to the previous result, 

(l + i;)Ve-2-t , 



(87) 



such that, for t +00, p diverges along every single 
u ray as e"'"^''*, because of the blueshift factor (1 -I- v)'^. 
Hence, p will diverge on the whole line , which is 
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a Cauchy horizon for the geometry (so this result does 
not contradict the FuUing-Sweeny-Wald theorem p^). 
We then deduce that the warp-drive spacetime is again 
hkely to become unstable in a very short time. 

As a closing remark, it is perhaps important to stress 
that this divergence of the RSET on the Cauchy hori- 
zon is of different nature with respect to the one found 
in Sec. IIV Bl In fact, the divergence of Sec. IIVBI is in- 
trinsically due to the inevitable transient disturbances 
produced by the formation of the horizon. In this sense 
it is a new and very effective instability present every 
time a white horizon is formed in some dynamical way. 
On the contrary, the just found divergence on can 
be seen as due to the well-known infinite blueshift suf- 
fered by light rays as they approach a Cauchy horizon, in 
this specific case as due to the accumulation of Hawking 
radiation produced by the black horizon. In this sense it 
is not very different from the often claimed instability of 
inner horizons in Kerr-Newman black holes JL6, 17, lj|. 

Of course the divergence and the appearance of the 
horizon would be avoided if the superluminal warp 
drive were sustained just for a finite amount of time. In 
that case, no Cauchy horizons would arise and no actual 
infinity would be reached by the RSET. However, the 
latter would still become huge in a very short time, in- 
creasing exponentially on a time scale 1/k w A/c, where 
A is the thickness of the warp-drive bubble. Note that, 
in order to get a time scale of even 1 s, one would need 
A « 3 X 10* m. 



V. SUMMARY AND DISCUSSION 

We have described the causal structure of both eter- 
nal and dynamical warp drives. In a geometric optics 
approximation we have studied the propagation of light 
rays in dynamical geometries and found that the same ex- 
ponential relation between affine coordinates on J'~ and 
is recovered at late times (large u) as in the case 
of black hole spacetimes. Given this relation it is un- 
avoidable the conclusion that indeed a Hawking flux will 
be observed by any observer inside the warp-drive bub- 
ble far from the black horizon. Indeed, the calculation 
shows explicitly the onset of such a flux. This radiation 
is produced at the black horizon and soon fills the in- 
terior of the bubble, traveling rightward at the speed of 
light. The central region of the warp drive behaves like 
the asymptotic region of a black hole: In both these re- 
gions the static contribution (pst) to the energy density 
vanishes so that the total energy density is due solely to 
the Hawking radiation generated at the black horizon. 

When creating a warp drive one not only forms a black 
horizon but also a white one. Both arc sources of right- 
going radiation. To understand better the nature of this 
radiation we have calculated the RSET in this geome- 
try and, in particular, the energy density as measured 
by freely falling observers. In this way we recover that 
the RSET docs not diverge at the horizons at any finite 



time. The singular behavior of the static terms (or vac- 
uum polarization terms) in the RSET at the horizons is 
canceled by the leading contributions of the dynamical 
right-going terms, or what is equivalent, by the presence 
of Hawking radiation at both the horizons. 

It is however easy to see that the subleading terms of 
Eq. ([551) behavior is rather different between the black 
and the white horizons. The subleading term in the 
RSET associated with the formation of the black hori- 
zon does not produce any significative backreaction on 
the horizon itself. In fact, this term is just a transient 
which decays exponentially.^ 

The formation of a white horizon is also associated 
with a similar subleading term, but this time it accu- 
mulates onto the white horizon itself. This causes the 
energy density p seen by a free-falling observer to grow 
unboundedly with time on this horizon. The semiclassi- 
cal backreaction of the RSET will make the superluminal 
warp drive become rapidly unstable, in a time scale of 
the order of 1/^2, the inverse of the surface gravity of 
the white horizon. Indeed, if one trusted the QI [alllj 
the wall thickness for vq k, c would be A < 10^ ip, and 
the surface gravity n > 10~^ tp~^, where tp is the Planck 
time.^" This means that the time scale over which the 
backreaction of the RSET would become important is 
T ^/^ ^ 10^ tp. Indeed, even forgetting about the QI, 
in order to get even a time scale r ~ 1 s for the growing 
rate of the RSET, one would need a wall as large as 3 x 10* 
m. Therefore, most probably one would be able to main- 
tain a superluminal speed for just a very short interval of 
time. In addition to the above mentioned growing term 
on J>^2~ '^c have shown that there is also an unbounded 
accumulation of Hawking radiation on ■ Also this 
contribution will very rapidly lead to a significant back- 
reaction on the superluminal warp drive and to some sort 
of semiclassical instability of the solution (that will most 
probably prevent the formation of the Cauchy horizons 
at late times). 

Interestingly, recent investigations (loj seem to im- 
ply that the above found asymptotic divergences of the 
RSET on the white horizon might disappear if Lorcntz 
symmetry gets broken at high energies. While in this 
case we still expect that the RSET will acquire large val- 
ues soon after the white horizon is formed, we do not 
know if this would be enough to prevent the sustainabil- 



^ Furthermore, the analogy with a black hole originated by a star 
collapse of [To| also allows us to infer that the warp drive must 
be created very rapidly in order to avoid a huge accumulation of 
vacuum polarization at the horizon and so a huge initial value of 
the energy density. 

" In 1 -I- 1 dimensions the warp-drive configuration is actually a 
vacuum solution of Einstein equations. In this case, QI will not 
impose any conditions on the size of the wall thickness. However, 
as we will show, we expect our results to be valid also in 3 -I- 1 
dimensions. Therefore, we use the wall thickness of 3 -I- 1 warp 
drives to obtain a realistic estimation for the Hawking tempera- 
ture and the time scale of the exponential growing of the RSET. 
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ity of the warp drive. We think that this issue could be 
subject for future research perhaps within the context of 
anafog models of gravity (see [20| for a complete review 
on analog models) where it might even be addressed ex- 
perimentally. 

Even if the above described semiclassical instability 
could be avoided by some external action on the warp- 
drive bubble (or by some appropriate UV completion of 
the quantum field theory, like in [l^), the QI lead to 
the conclusion that the Hawking radiation in the center 
of the bubble will burn the internal observer with an 
excruciating temperature of Th ~ k > 10"^ Tp, where 
Tp is the Planck temperature, about 10'^^ K. This would 
prevent the use of a supcrluminal warp drive for any kind 
of practical purpose. If we do not trust the QI, this high 
temperature might be avoided by making thicker walls. 
For instance, with A ~ 1 m, one obtains a temperature 
of about 0.003 K (roughly the temperature of radiation 
at a wavelength of 1 m).^^ 

Finally we want to comment on the fact that in this 
paper a l-fl calculation was performed. Generally in 
spherically symmetric spacetimes this could be seen as 
an s-wave approximation to the correct results. How- 
ever, this is not the case for the axisymmetric warp-drive 
geometry. However, we do expect that the salient fea- 
tures of our results would be maintained in a full 3+1 
calculation (most probably a numerical one) given that 
they will still be valid in a suitable open set of the hori- 
zons centered around the axis aligned with the direction 
of motion. 

In conclusion, we think that this work is casting strong 
doubts about the semiclassical stability of supcrluminal 
warp drives. Of course, all the aforementioned problems 
disappear when the bubble remains subluminal. In that 
case no horizons form, no Hawking radiation is created, 
and neither strong temperature nor white horizon insta- 
bility is found. The only remaining problem is that one 
would still need the presence of some amount of exotic 
matter to maintain the subluminal drive. 
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However, some very effective taming for the growing backreaction 
at the wfiite horizon would be needed also in this case, given the 
previously estimated RSET growing rate. 



APPENDIX A: AN EXAMPLE OF DYNAMICAL 
WARPDRIVE 

In this Appendix we want to present a model of dy- 
namical warpdrive, which we used for numerical calcu- 
lation. It is particularly useful because it allows one to 
carry out an almost complete analytical treatment. We 
adapt here the method presented in [13, [2l| for stellar 
collapses to the creation of a warp drive from Minkowski 
spacetime. We can choose a simplified piecewise velocity 
profile, which has the relevant properties of a dynamical 
warp drive, i.e., it describes a flat geometry at early times 
and coincides with the metric in Eq. ^ after some finite 
time t > 0. 

Using our velocity profile defined in Eq. ([5]) 



v{r) = ac[/(r) - 1] , 



(Al) 



we can define a dynamical profile by replacing ij in 
Eq. (inil) by Wkink 



t'kink(f, = 



v{r) 



if 
if 



\r\>m 



(A2) 



where S,{t) is a monotonically increasing function of t, 
such that £,{t) 0, for t — oo, and ^(^h) = ^2 = — ri. 

One may wonder whether defining a velocity profile 
with a kink, as in Eq. (jA2[) . may lead to unphysical phe- 
nomena. Indeed this computational trick induces some 
spurious effects, but these features arc just transients and 
do not affect the results at late times. 



1. Right-going rays 

We apply the same procedure of [l3] to calculate the 
exact relation between the past null coordinate U in 
and the future null coordinate relevant at ■ To find 
this relation one has to find the integral curves of the ray 
differential equation ([55)) , which becomes at early times 



dr 
'dt 



and in a neighborhood of Jfft {t +oo and 
Fig. [3]) 



(A3) 



r2, see 



dr 
'dt 



+ v{r) = K{r2-r) + ((ra - rf) . (A4) 



Integrating the first equation we obtain the obvious result 



t = C 



and for the second one 



t = D In (r2 - r) . 



(A5) 



(A6) 



Following [TJl we identify initial events P = {ri,ti), with 
fi ~ cti, final events Q = {rf,tf) with rj ^ 7*2 — e""*^, 
and intermediate events O = (ro,to). 
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Let us define 



Using tiie definition of U in Eq. (jA7[) we obtain 



U 



lim U 



lim 

t f — > + oo 



tf + -\n (ra-r/) 



(A7) 
(A8) 



Integrating Eq. (j33p between P and an intermediate 
event O in the in region of this spacetime (where the 
velocity profile depends only on i), we find 



f - n 



dtic+mm 



(A9) 



c c 



to 



, (AlO) 



Then, integrating between an intermediate event O 
now in the out region (where the velocity profile depends 
only on r) and Q we find 



tf -to - 



'■^ dr 



(All) 



and using the definition in Eq. (jM]) 



U = tn 



lim 



- In (r2 - Tf) + 

K 



dr 



c + v(r) 



to + - In [r2 - ro] - - In [ro - ri] + - In [r-z - ri] + hm 

K K K 



dr 



1 



1 



1 



c + w(r) K (r2 — r) K(r — ri) 



(A12) 



It is easy to see that the limit in the previous expression is finite for whatever ri < ro < r2 so that wc end up with 
the relation 



1 1 1 r 

w = to + - In V2 -ro\ In [ro - ri] + - In [r2 - r^\ + / 

K K K J J. 



dr 



1 



1 



c + w(r) K (r2 — r) K(r — ri) 



(A13) 



We now want to find the relation between U and u. It 
is possible to find a particular form for the kink function 
^ (see Fig.|4]and its caption for an example), such that all 
the rays cross either the left kink r = —£,{t) or the right 
kink r = +^(t) only once. The event "crossing the kink" 
can be seen as belonging to both the in region and the 
out region. Therefore the relation U = p{u) can be found 
by eliminating Iq and rg between expressions (jAlOp and 
(|A13p taking into account that ro will be either — ^(to) 
or +£,{to)- Figure m shows an exact numerical calculation 
of ray propagation, performed with this method. 



However, in order to study the late-time behavior of 
both particle production and the RSET, we only need the 
previous relation (see Sec. IIIip for large and positive val- 
ues of u (the location of the black horizon) and for large 
and negative values of u (the location of the white hori- 
zon) . We want to show here that our particular example 
leads in fact to the general relation Eq. ([40|) . In this spe- 
cific case, taking the limit for |u| — > 00 corresponds to 
study the previous relation for light rays crossing the left 
kink at ro = —£,{to), when ro is very close to ri [respec- 
tively, crossing the right kink at ro = +^(to)j when ro 
is very close to r2]. First note that, when t ^ tn (the 
time of the first appearance of the trapping horizons) the 



function ^ can be expanded in the following form: 

m=r2 + \{t-tH) + o[{t-tHf] , (A14) 

-m = ri-\{t-tH) + o[{t-tHf'\ , (A15) 

with A a positive constant (remember that we have cho- 
sen for simplicity ri — — r2). Defining 



C^BH — iff ± 

WH 



it is easy to see that 

C/± = C/bh T-{tH-to) 

WH C 



O 



dtvim) 



(A16) 



(A17) 



where with C/-|_ wc denote the form of the function U{to) 
for rays crossing the kink close to ri [respectively, with 
U- we denote the form of the function C/(to) for rays 
crossing the kink close to to ri]. By expanding Eq. (jA13P 
in the same limit to — > tn and retaining only the domi- 
nant term we obtain 



u± ~ iF-ln[A(tH - to)] 



(A18) 



where we define u± in the same fashion as U±. Putting 
together the last two results we obtain (inside the bubble) 
the general relation Eq. (|40|) : 

U(u ±00) = Ubh T Aie^F-^". (A19) 

WH 
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2. Left-going rays 

In Sec. IIVI we stated that it is possible to find a wide 
class of transitions from a Minkowskian geometry to a 
warp-drive one such that the contributions to the RSET 
due to left-going modes are just transient, vanishing at 
late times. We show here that this is the case for our 
specific model. 

In order to do so. we shall need the relation between 
ingoing and outgoing left-going rays. Left-going rays are 
the solution of the differential equation (pij) which be- 
comes at early times 



dr 
'dt 



(A20) 



while inside the bubble we cannot take any limit as in 
Eq. (|A4[) because left-going rays are not confined inside 
the bubble, but escape through the black horizon, as we 
already noticed in Sec. IIIIl Integrating the equation at 
early times wc obtain 



t = C 



(A21) 



We identify initial events P = (ri,ti), with cti, final 
events Q = {rf, tf), and intermediate events O = (ro, to) 
either in the in region (when the velocity profile depends 
only on t) or in the out region (when the velocity profile 
depends only on r). 

Using w defined as in Eq. (|^^ . choosing arbitrarily the 
constant of integration, we obtain 



u!{tQ,ro) = to 



dr 



v{r) 



(A22) 



For W we proceed as for U. We integrate Eq. ([H 
between P and an event O in the in region 



to 

ro -r, = I dt[-c + v{at))] 



Defining W as 



W = hm {U + - 

ti—t — ca \ c 



(A23) 



(A24) 



we obtain 



ro 1 



to 



dtvm) ■ (A25) 



W = q{w) can be found eliminating 
and (jA25p in the limit where {to, ro) 



Now, the relation 
to, ro from Eqs. (H^ 
is on the kink. 

We can now show that indeed the term pdyn-to defined 
in Eq. (|74p is just a transient for this model. We see from 
Fig. [4] that, as time grows, points inside the bubble are 
reached by w rays that intersect the kink later and later. 
Using Eqs. gH) and (|M5l) we can estimate q, calculating 



the derivatives of w and W with respect to t on the kink, 
at crossing time t^- 



dW 


dW 


1 dw 


dw 


~ ~di 




(■■ 


i{tc) 


mtc))^ 1 




c 








_ ^ mtc)) 



c-mtc))^ 

v{+oo) 



1 - 



(A26) 



This is a constant value greater than 0. As a conse- 
quence, all the derivatives of q go to zero as time grows 
and pdyn-w must go to zero in the same way. In con- 
clusion, the dynamic term originated by the distortion of 
left-going rays can be different from zero when the hori- 
zon is created. However, this is only a transient term that 
is brought toward , i.e., outside the bubble. That is, 
it is possible to create a transition region such that some 
radiation traveling leftward is produced only at the onset 
of the warp drive, but there is not any phenomenon like 
Hawking radiation originated in this way. 



APPENDIX B: NULL COORDINATES 

In this Appendix wc summarize some results concern- 
ing systems of null coordinates. In Sec. lIIII we defined two 
sets of null coordinates {U, W) and {u, w) obtained by in- 
tegrating equations of light propagation, respectively, in 
the in region, in which the metric is Minkowskian, and in 
the out region, in which it depends on r (after the kink, 
in the particular example of Appendix [S| . Now, let us 
rewrite the metric using the set of coordinates {u, w). In 
the out region, using Eqs. ([M)) and 



1 



c + v{r) ' 



c — v{r) 



(Bl) 



Ut 



Wt 



Now we 



put these expressions into Eq. (j45|) : 



ds^ = —C{u,w)dudw 

~ -~C{u, w) {utdt + Urdr) (wtdt + Wrdr) 
C{u, w) 



'^%r\c^dt^ ~[dr~v{r)dtf\ . (B2) 



;(r) 

By comparison with Eq. (j3ip we obtain 
C{u, w) ~ (? ~ v(r)^ , 



(B3) 



Note that the definition of U and W can be extended 
in the out region by following the light rays coming from 
the in region (in the model of Appendix [X] by matching 
on the kink the solutions of U and W with those of u 
and w). For instance, one can naturally define U{t,r) in 
the out region as U{t,r) = p(u(t,r)) and in analogous 
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ways for the other null coordinate. This means that it 
makes sense to take derivatives of the null coordinates 
with respect to r and t in all the spacetime. Wc write 
here the derivatives of U and W in the out region: 



Ur = p(u)Ur 



p{u) 



c + v(r) 
Ut =p{u)ut ^p{u) , 



Wr = q{w)Wr 



c — v{r) 
Wt = qiw)wt = q{w) . 



(B4) 



APPENDIX C: ASYMPTOTIC EXPANSION OF 

/(«) 



where T'± are analytic functions. We now expand it up 
to the third order: 



U =p(u^ ±00) = Ubh + Aiie^""" 

WH 

A 



A2± 



where the upper (lower) signs correspond to u ^ +00 
{u —>■ —00). For simplicity we have defined the coeffi- 
cients Ai+ = —A^ and Ai- = A-. In order to calculate 
the stress-energy tensor we need to calculate the term 
Eq. ^ 



We know from Sec. IIIII that U can be expanded for 
large absolute values of u in a Taylor series in e^''". In 
fact we found that 



/(") 



3p2 - 2p P 



(C3) 



U± = p{u ±00) = P±(eT'^") 



(CI) 



We have 



p{u) = K2^i±e=F™ 



Ai± 
2A2± 

4A 



2A 



i± 



2± 



Ai-_ 



2Ai± 



so that 



(C4) 
(C5) 
(C6) 



/(") 



K4A2^eT2«« |l + 2(A2±Mi±)eT'=" + 


4(^2±Ml±)' 


- ^3±Ml± 






K2^2^eT2«« |l + 2(A2±/Ai±)eT™ -|- 


"(A2±Mi±)' 


+ ^3±Ml± 




-0(e=F3™)| 



^2 < 1 



4(A2±Mi±)'-A3±Mi± 



1 + 2(A2±/Ai±)eT« 



gT2K« ^_ Q (e=F3«") I X <^ 1 



(A2±Mi±)'-f A3±Mi± 



k2<^1 



1 + 2(A2±Mi±)eT« 
2 



e=F2«"-(-C'(e 



\Ai±J Ai± 



gT2«« ^ Q (e=F3««) ^ , (C7) 



and we can finally use the expansion of Eq. (j36p . valid for points close to the horizons, to obtain 



fiu) + 



^2± 

Ai± 



A 



3± 



Au 



eT2Kt (r ~ n.2f + O [{r - n.2f) 



(C8) 
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